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TJ} ' Abstract 

, Sine-Gordon model with variable mass (VMSG) appears in many physical systems, ranging from 

■ the current through nonuniform Josephson junction to DNA-promoter dynamics. Such models are 

usually nonintegrable with solutions found numerically or peturbatively. We construct a class of 
Q\ ■ VMSG models, integrable both at classical and quantum level with exact soliton solutions, which 

t-H , can accelerate, change their shape, width and amplitude simulating realistic inhomogeneous systems 

at certain limits. 

m ■ 

PACS: 05.45.Yv,03.65.Fd,11.10.Lm,11.55.Ds 

Sine-Gordon (SG) model enjoys a special status among nonlinear integrable systems for its inherent 
richness and wide range of applications in different fields [1, 2, 3, 4, 5, 6, 7, 8, 9]. Apart from the 
fascinating properties of an integrable system, the SG model exhibits relativistic invariance and integer- 
valued topological charge represented by solutions like kink, antikink, breather etc. [10], together with 
the quantum integrability described by the Yang-Baxter equation (YBE), which for the SG leads to 
, quantum algebra su q {2) [11, 12]. 

Solitons in the standard SG model, as in other integrable systems, move with constant velocity and 
shape. In realistic situations however due to inhomogeneity of the media solitons may exhibit more 



complex motion with changing velocity and shape [4, 6, 7], which may be used also as a desirable 
effect for fast transport, fast communication, or even for a possible soliton gun [8]. In particular 
inhomogeneity can lead to SG models with variable mass (VMSG) in describing fluxon or semi-fluxon 
dynamics in Josephson junction (JJ) with impurity or nonuniform critical current [2, 6], spin wave 
propagation with variable interaction strength [5] , DNA-promoter dynamics in nonuniform background 
etc. However such inhomogeneities usually destroy the most cherishable property of the SG model, 
i.e. its integrability and hence the solutions can be extracted only numerically or at best peturbatively 
[2,4,6,7]. 

We observe that, though the integrability of SG model is spoiled by a variable mass, or a variable 
velocity, it can be restored if both of them vary simultaneously following certain rule. Therefore we 
can construct a VMSG model, integrable at the classical and the quantum level, allowing analytic 
soliton solutions. Such exact solitons nevertheless show intriguing accelerated motion with changing 
shape, amplitude and width, simulating thus realistic inhomogeneous systems [4, 6, 7] and describing 
them analytically at certain limits (Fig. lb-d). 

To clarify our strategy we focus on the linear spectral problem of the SG model: & x (x,\) = 
U(\,x)${x,\), 

&t(x, A) = V(X, x)Q(x, A), with its Lax pair [13]: U = | (— uta 3 + mk\ cos ^cj 2 — mho sin fc 1 ) , and 



1 



V = I (— u x a 3 — mko cos ^a 2 + mk\ sin ^a 1 ) , where &o(A) = 2A + j^, k\(X) = 2A — with spectral 
parameter A. Compatibility condition (E^t = <&tx or the related flatness condition leads to the SG 
equation, for constant mass mo and spectral parameter Ao- Recall that in the inverse scattering (IS) 
method solitons are obtained at discrete spectrum A n , n = 1,2, ... N, with velocities of SG solitons 
linked to these values of the spectral parameter. Therefore variable soliton velocity should have a 
variable spectral parameter A, which alone however violates the compatibility condition. Fortunately 
by making mass m also a space-time dependent variable, we can get the VMSG equation 

u tt ~ u xx + m 2 (x,t)smu = 0, , (1) 

with the constraint: {kom)t + (mki) x = 0, {kim)t + (mko) x = 0, which can be reduced to two free 
field equations: 

Ktt ~ Kxx = 0, ptt - Pxx = 0, for k = Inm(x.t), p = lnX(x,t). (2) 

Note that the set of equations (1-2) is a new integrable relativistic system, generalizing SG equation, 
and a reduction (at the free field limit of spectral dilatation field p) of the conformal affine Toda model 
(CATM) [14]. However since we are interested here in application to physically relevant inhomogeneous 
models, we would consider k, p as given inhomogeneous functions by restricting to particular solutions 
for variable mass and spectral parameter: 

m(x,t) = m /+/-, \(x,t) = ^oj~, (3) 

with f± arbitrary smooth functions of x± = x ± t, respectively. Thus we obtain an integrable SG 
equation (1) with variable mass m(x,t) = rriof+f-. Note that due to explicit space-time dependent 
coefficient, it is no longer a relativistic and translational invariant model. However if we demand such 
invariance, we simply get back the constant mass SG model, as shown in [15]. We can control mass 
m(x,t) in the VMSG model (1) by choosing suitably the inhomogeneity functions f± for different 
physical situations, showing a variety of soliton dynamics as in Fig. 1 a-d. 

For obtaining exact solutions of the VMSG model (1), we can apply well known methods de- 
signed for integrable systems [10], e.g. Hirota's bilinearization and the IS formalism. Hirota's 
method for soliton solution for the standard SG equation is given through a cleverly chosen ansatz 
u = — 2iln with g± as conjugate functions, which converts the SG equation into a bilinear form, 
admiting solution for g± as expansion in plane-waves. For SG model (1) with variable mass and 
velocity the same ansatz seems to work, only the plane waves should be replaced by their generalized 
form: = ^. e U x (^,x,t)-T(X n ,x,t)) j where X (X n ,x,t) = J x dx'm{x' ,t)k 1 {\ n ,x' ,t), T{X n ,x,t) = 

ydt'm{x,t')k {\ n ,x,t'). This gives the exact soliton solutions of (1) through the expansion: g± = 
l±gW, for 1-kink and g± = 1 ± (gW + g^) + s(i(0i - # 2 ))5 (1 V 2) , A a = §e e % a = 1, 2 for 
2-kink etc., with the scattering amplitude s(9) = tanh 2 9, while for A2 = — AJ = r]e 10 , one gets the 
kink-antikink bound state (breather solution). 

Similarly we can apply the IS formalism [10] to (1), for which the crucial step is to use the analyticity 
of Jost solutions <I>, based on their behavior at A — > ±00. This holds equally for the inhomogeneous 
extension of the SG model, where the asymptotic plane waves should again be replaced by their 
generalized form. 1-kink soliton with Ai = if] can be obtained explicitly, either from the Hirota's or 
from the IS method as 

u = 4tan~ V), C = \{X(i-q, x, t) - T{i V , x, t)), (4) 

with variable soliton velocity v s (x,t) = = klfyx'f) • Kink solution (4) gives a localized soliton for 
sin I = co l h( . , which are actually shown in Fig 1. 
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To see the effect of various inhomogeneities on the soliton dynamics, we consider some concrete 
integrable cases. Notice that variable mass i) mo(x 2 — t 2 ) n , invariant under relativistic motion, yields 
( for n = 1) the exact kink solution (4) with ( = 1 j(2r](x — i) 3 + ^-(x + t) 3 ), the evolution of the 
corresponding soliton is depicted in Fig la. The intriguing change in the soliton shape, width and 
velocity during its motion is clearly seen. Position-dependent mass can be achieved in this case at 
t —> and therefore phenomenon like fluxon propagation through JJ with local defect mox 2 may be 
described with the above analytic soliton solution at a short time interval limit, as shown in Fig lb. 

Other forms of integrable VMSG equations can be obtained for mass ii) y/2mo cos a q(x ± t), a 
being an arbitrary parameter. We derive for the first case (with a = 1) kink solution (4) with 
C = mo(ko(r])x — k\{rj)t + ^ sin 2q(x + t)), having soliton velocity v s = mod{k\{n) — ^ cos 2q(x + t)) 
and width d = (mo(fco(?j) + ^ cos2q(x + t))) -1 , both of which oscillate periodically in space-time, as 
evident from Fig. lc. Notice that variable mass of this type is particularly interesting, since it can 
describe an important physical situation, namely the SG equation parametrically driven by a plane 
wave [17]. 

One can get a similar integrable case with mass iii) mo{2 cos q{x + t) cos q(x — t)) ~, which at short 
time interval limit (t — > 0) gives ~ m{x) = tHq (cos qx) a , while for evolution limited to a small space 
interval (x — > 0): « m(t) = rho(cosqt) a . Recall that a physically motivated spin chain with coupling 
constant changing periodically in space, can be described by a VMSG with mass m(x) = mo(cos qx) a , 
where a = 2^k> with K > \ being an important parameter of the system [5]. Similarly a real 
oscillator chain pumped by an alternating current [16] can be linked to a VMSG with mass m{t) = 
(cos qt) 2. Therefore we may conclude that the analytic solution of our VMSG equation can describe 
the inhomogeneous spin wave dynamics [5] or evolution of forced oscillators [16], at least at short 
time or space interval limit. Alternatively, this realistic spin (or oscillator) model can be tuned to the 
integrable VMSG with mass mo cos" q(x + 1), by making its coupling strength oscillate periodically 
also in time (or space). 

In most physical situations the inhomogeneity of the media leads to VMSG equation, with only 
position-dependent mass m(x). Therefore we explore to find, when such equations can be integrable 
in the entire space-time and conclude from our result (3), that the VMSG becomes integrable only 
for the space-dependent mass iv) m(x) = moe p ( x ~ x °) with p = const., which explains also why most 
of the realistic VMSG models having a different position-dependent mass (e.g. [5]) turn out to be 
nonintegrable systems. Exact kink solution for the integrable m(x) is obtained from (4) as u = 
4tan _1 (e 1 '), ( = -ko(t)m(x), m(x) = exp(/o(x — xq)), ko(t) = cosh(# — p(t — to)). The corresponding 
soliton velocity and width are v s = tanh(9 — p(t — to)) , and d = (m{x)ko{t))~ l , showing how the soliton 
shape changes and how it accelerates, decelerates or can exhibit boomeron [18] like property. This 
scenario is close to the predicted behavior of solitons in the dynamically active promoter zone of the 
T7Ai DNA [4]. Notice that for p > 0, since m(oo) = oo, m(— oo) = 0, the kink solution yields u(oo) = 
2tt, u(— oo) = 7r , and hence corresponds to the topological charge Q = ^(u(co) — u(— oo)) = ^. This 
intriguing fact might serve as an integrable theory for the semi-fluxon, observed in unconventional JJ 
[19]- 

At p — > 0: ( — > Co = m o(ko (x — xq) — k\ (t — to)) and the standard SG soliton is recovered. 
Therefore we can acess the solitonic behavior in realistic models with any mass deviation from its 
constant value, with high degree of accuracy, by approximating through expansion in powers of p. 
Fig. Id shows that a static soliton remains static, when placed in a region with constant mass while 
an initially static soliton can move with accelerated (or decelerated) motion, when placed in a zone 
with variable mass, resembling closely the scenario of the VMSG soliton in the DNA chain, which with 
zero initial velocity in an inactive regions (with constant mass due to almost uniform background of 
two types of base pairs) remains static, while in the active promoter region with variable mass (due 
to significant difference in the number of lighter (A-T) and heavier (G-C) base pairs) the same static 
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Figure 1: Exact soliton solutions: s = sin t) of the integrable VMSG equation with variable mass 
m. a) m = mo(x 2 — i 2 ) having an intriguing flattening of the soliton at the center, b) Short time 
interval limit of the above soliton showing the flattening prominently, c) m = 2m.Q cos q(x + t) with 
oscillatory behavior of the soliton. d) Static soliton in the zone [x < 1.2) with m = const and initially 
static soliton in the zone (at x = 4.8) with variable mass: m = rriQ exp(px) ~ m(l + px), with p = 0.1 
moves backward with acceleration, resembling soliton propagation in inactive/active promoter region 
in a DNA chain. 

soliton can acquire rich accelerated motion [4]. 

Finally we intend to show that the integrable VMSG model constructed here can be raised 
to the quantum level and the algebraic Bethe ansatz (ABA) developed for the constant mass SG 
model [11] can be adopted successfully for it. Quantum lattice regularized SG matrix Lax opera- 
tor Uj(X,Sj,m), j = 1,2, ...,L involves quantum-spin operators Sj(uj), (uj , pj , m) , expressed in 
canonical operators Uj, Pj = Uj and the mass parameter m, which should be generalized now to site 
dependent parameter rrij [20]. Note that, the trigonometric i?(^)-matrix associated with our quantum 
integrable VMSG model remains the same, since it depends on the ratio of two spectral parameters, 
in which x, t-dependence enters, as seen from (3), only multiplicatively and therefore gets canceled. 
Moreover, YBE being a local algebra (at each lattice site j) is not affected by inhomogeneity and 
yields the same quantum algebra su q (2); only with the replacement of constant m by a site-dependent 

function rrij m its structure constant: [bj, b k \ = Oj^rrij s i na • 

The aim of ABA is to solve exactly the eigenvalue problem of trT(X) = A(X) + D(X), with 
T(A) = Uj Uj(X), generating all conserved operators including the Hamiltonian, with the eigenstates: 
|Ai,...,A n >= YYJ B(Xj)\0 >. Ti2(A) = B(X) acts as creation operator, while T2i(A) = C(A) as 
destruction operator annihilating the pseudovaccum: C(A)|0 >= 0. Following closely [11], but gener- 

(2) 

alizing for the site-dependent mass rrij, we can construct the local pseudovaccum |0 >= JTj \ ftj >, a 
crucial step in ABA, by combining the action of consecutive pair of Lax operators: UjUj+i\Q > [21]. 
Consequently the vacuum eigenvalues are generalized for the quantum VMSG model as j4(A)|0 >= 
a (m) |0 >, £>(A)|0 >= P (m) \0 >, where a {m) = Uj a(9, ^), /3 (m) = Uj a* (6, ^) with a(9, ^) = 
^ + <5 2 mim2(cosh(2# + ia)). This yields the exact eigenvalue for the conserved quantities: trT(X) 
as A(A; Ai, . . . , An) = a {m) f|" f(-f) + P( m ) Uj /(^)) where /(^) is expressed through the ele- 
ments of the -R(^)-matrix, which remains unchanged for the VMSG model. The Bethe equations for 
determining the parameters {Xj} are generalized similarly by taking m — > rrij. 

Since our VMSG model can be reduced from CATM, a coordinate transformation (x, t) — > (X, T) 
exists, which can take the SG model with variable mass formally to the SG model with constant 
mass, though the domain might shift to an unphysical region and singularities might arise. Such a 
nonlinear transformation, amounting to going to a noninertial frame, takes particularly simple form 
in the light-cone coordinates : X± = J dx±f± [15, 14]. However for investigating the physical effect of 
a given inhomogeneous medium inducing accelerated and shape changing solitons, one has to analyze 
the model in its original form with variable mass. Similar situation arises also in other inhomogeneous 
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systems with integrable nonisospectral flow, e.g. [22] in the study of accelerated solitons in plasma 
through NLS equation, in discrete NLS, in Toda chain, in matrix Schrodinger problem with boomeron 
solution [18] etc. In most of these models though the inhomogeneities could be removed by tricky 
nonlinear transformations, the investigations were carried out in the original systems due to their 
physical relevance. Surprisingly this long list of literature devoted to various inhomogeneous integrable 
models does not include the well known SG model and also avoids any quantum treatment, except 
perhaps a recent attempt [23]. This enhances therefore the importance of the present result, which 
explores the inhomogeneous integrable SG model, presents its exact solution both in the classical 
and in the quantum case. Its analytic soliton solutions can simulate at certain limits realistic events 
like, fluxons in nonuniform Josephson junction, dynamics of spin waves with variable coupling, DNA 
solitons in the active promoter region etc. Regulating the position-dependent mass in integrable 
VMSG model one can create semi-kink solution suggesting a possible exact theory for semi-fluxon. 
The author likes to thank Profs. P. Mitra (SINP) and D. Sen (IISC) for helpful discussions. 
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